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INTRODUCTION 

This  manual  has  been  prepared  to  aid  field  personnel  in  using  the  Bureau 
of  Land  Management's  (BLM)  Geographic  Information  System  (GIS)  to  store, 
retrieve,  analyze  and  model  spatial  data  for  use  in  coal  tract  analysis. 
GIS  can  be  used  to  generate  contour  maps,  calculate  coal  tonnage  and 
overlay  several  maps  to  create  a  composite  map.  GIS  can  also  be  used  to 
edit  data  files  and  generate  reports. 

The  two  GIS  software  packages  discussed  in  this  manual  include  the  Map 
Overlay  Statistical  System  (MOSS)  and  the  Map  Analysis  and  Processing 
System  (MAPS)  . 

If  you  have  any  GIS-related  questions,  contact  the  HOTLINE  at  FTS  776- 
0990  or  commercial  (303)  236-0990. 

SECTION  ONE:  METHODS  TO  CREATE  A  DRILLHOLE-LOCATION  MAP 

GIS  can  process  non-map  data  and  transform  the  data  into  a  positionally 
correct,  coordinate-referenced  map.   To  create  the  drillhole-location 
map,  there  are  four  general  methods  in  which  this   an  be  done: 

Method  1:  DIGITIZING 

A  drillhole-location  map  can  be  generated  by  digitizing  points  from  a 
basemap  using  an  electronic  graphics  tablet  and  a  manually-operated 
cursor  that  encodes  the  coordinates  for  each  point  using  the  Automated 
Digitizing  System  (ADS) .   The  digital  coordinates  can  then  be 
auromatically  reformated  into  a  MOSS  map  that  is  compatible  with  multiple 
attributes.   (Scan  digitizing  may  be  economical  in  the  near  future.) 

Method  2 :  COALSYS 

If  the  data  are  in  COALSYS,  use  that  software's  programs  to  generare  an 
import  file  for  MOSS.   The  COALSYS  manual  will  explain  how  to  do  this. 
This  method  generates  an  import  file,  not  a  map.   The  import  file  is  then 
used  with  the  IMPORT  command  to  generate  a  map. 


Method  3:  GGWL 

A  program  called  Generate  Graphic  Well  Locations  (GGWL^  can  generate  the 
UTM  coordinate  values  from  location  footages  (feet  fror..  section  lines)  . 

Method  4;  XYSUBJECT 

I,f  the  data  are  from  some  other  data  base,  download  an  ASCII 
(alphanumeric  text)  file  that  contains  the  drill-hole  identifier  and 
coordinates  for  each  point.    Then  use  the  MOSS  UTILITY-XYSUBJECT  command 
(menu  option  #16)  to  reformat  the  coordinate-data  file  into  an  import 
file  that  can  be  used  by  the  IMPORT  command  to  generate  a  map. 

Method  4  (XYSUBJECT)  is  addressed  in  the  following  sections  of  this 
manual . 

SECTION  TWO:  USING  XYSUBJECT  METHOD  TO  CREATE  THE  IMPORT  FILE 

This  is  a  short  example  of  a  coordinate  file: 

DH-1    44.6789  107.7654 
DH-2    43.9999  107.5555 

XYSUBJECT  asks  if  the  coordinates  are  in  decimal  degrees,  meters  or  feet 
and  displays  a  counter  line  above  the  first  line  in  the  coordinate  file 
as  shown  below: 

1234567890123456789012  3  4567899123  45678901234567890 
DH-1    44.6789  107.7654 

The  program  then  asks  the  user  enter  the  "format  statement".   The 
provided  format  is  unique  to  the  coordinate  file.   In  this  case,  the 
format  for  the  above  example  is: 

(T16,F8.4,T8,F7.4,T1,A4) 

This  instructs  the  program  to: 

tab  over  to  column  16  for  the  X-coordinate,  which  is  a  real  number 

with  a  field  length  of  8  spaces  and  4  spaces  to  the  right  of 

the  decimal  point 
tab  over  to  column  8  for  the  Y-coordinate,  which  is  a  real  number 

wit   a  field  length  of  7  spaces  and  4  spaces  to  the  ri -'h :.  of 

the  decimal  point 
tab  over  to  column  1  to  find  the  subject  which  is  character  in  cype 

and  has  a  field  length  of  4  spaces 

The  import  files  should  be  named  with  the  suffix  .IMP  so  that  they 
can  be  easily  identified  as  import  files.   XYSUBJECT  will  automatically 


name  the  output  file  XY.IMP.   This  should  be  renamed  as  the  next  run  of 
XYSUBJECT  will  overwrite  this  file. 

FACTS  TO  REMEMBER: 

The  subject  and  coordinates  in  the  coordinate  file  may  be  in  any 
columns  as  long  as  there  is  one  row  per  drill  hole. 

The  coordinates  may  be  in  latitude/longitude  (decimal  degrees) ,  feet 
or  meters. 

Latitude/ longitude  must  be  in  decimal  degrees.  Latitude/longitude 
coordinate  data  can  be  converted  to  decimal  degrees  using  the 
formula: 

decimal  degrees  =  degrees  +  (minutes  *  60)  +  (seconds  *  3  600) 

FIGURE  1  shows  an  example  of  an  output  file  from  the  XYSUBJECT 
command.   The  format  of  the  output  file  is: 


Line  1 ; 


15 


lOx 

3  0A 

5x 

15 


Item  number;  this  is  a  negative  integer  if 
coordinates  are  in  latitude  and  longitude; 
right  justified 

Ten  blank  spaces 

Subject,  e.g.,  the  drill-hole  number 

Five  blank  spaces 

Number  of  coordinate  pairs;  integer;  right 
justified 


Line  2  for  coordinates  other  than  latitude/longitude; 

FIX. 2    X-coordinate 
FIX. 2    Y-coordinate 

Line  2  for  latitude/longitude  coordinates: 

FX0.5    Longitude  in  decimal  degrees 
FX0.5    Latitude  in  decimal  degrees 


SECTION  THREE:  CnEATE  THE  DRILLHOLE-LOCATION  MAP 

STEP  1.  I^'"   AT  THE  THE  COORDINATE  T  .  IMP)  FILE 

FACTS  TO  REMEMBER: 

Input  file  must  be  in  MOSS  import  format. 

Must  enter  the  correct  number  of  map  subjects  that  are  in 

the  .IMP  file  to  be  imported.   To  find  this  number. 


FIGURE  1.   Example  of  an  import  file  that  include  identification 
labels  for  the  points  and  their  X-Y  coordinates. 


1  SJ45  1 
221484.00  4010740.00 

2  A-11  1 
210176.19  4015977.00 

3  B94  1 
210055.59  4014987.00 

4  A-12  1 
211221.59  4016059.00 

5  A-13  1 
211809.59  4015943.00 

6  A-14  1 
212393.38  4016064.00 

7  B69  1 
211936.88  4014941.00 

8  B50  1 
213124.88  4014932.00 

9  B39  1 
213855.78  4015471.00 

10  B23  1 
214769.59  4016115.00 

11  B25  1 
214716.50  4015469.00 

12  TL-7-1  1 
215810.69  4015678.00 

13  TL-7-2  1 
216508.50  4015585.00 

14  A-19  1 
215757.69  4014855.00 

15  A-20  1 
216330.78  4014789.00 

16  A-21  1 
217359.38  4014940.00 

17  TL-8-1  1 
217506.69  4015734.00 

18  TL-8-2  1 
218451.50  4015730.00 

19  A-23  1 
219150.09  4015748.00 

20  A-22  1 
218454.33  4014711.00 


either  SUST  the  .IMP  file  at  the  PRIMOS  level  or  use 
the  EMACS  editor  to  find  the  last  sequence  number. 

Must  enter  a  scale  factor:  latitude/longitude  files  will  use  100000 
as  a  scale  factor;  for  other  projections  use  100 „ 

Must  enter  the  correct  data  type  (point) . 

Other  information  asked  for  by  the  command  is  for  the 

map  header  information  and  has  no  effect  on  the 
calculations. 

FIGURE  2  is  an  example  of  the  IMPORT  command  in  which  the  import  file 
name  is  COAL. IMP  and  the  coordinates  are  in  latitude/longitude. 

STEP  2.  REPROJECT  THE  LATITUDE/LONGITUDE  MAP 

If  the  original  coordinate  file  was  in  latitude/longitude,  use  the 
PROJECTION  command  to  reproject  the  imported  map  to  another  projection, 
such  as  UTM  or  State  Plane.   FIGURE  3  gives  an  example  of  the 
PROJECTION  command  to  change  the  latitude/longitude  projection  into  UTM. 

SECTION  FOUR:  LINK  MULTIPLE  ATTRIBUTES  TO  THE  DRILLHOLE-LOCATION  MAP 

STEP  1.  CREATE  THE  DATA  FILE 

The  data  file  contains  the  multiple  attributes  (e.g.,  elevation)  for  each 
drillhole-location  on  the  map.   Naming  convention  gives  the  suffix  .DAT 
or  .ATT  for  the  file.   The  following  are  examples  of  multiple  attribute 
filenames: 


COAL. DAT 


MILS. DAT 


UINTAH. ATT 


PI. ATT 


The  multiple  attribute  f  i... .;  may  be  a  considered  a  two-dimensional  table 
composed  of  rows  and  columns.   Rows  in  the  table  represent  individual 
drillhole-locations  or  items: 

Row  1  will  correspond  to  Item  1  in  the  map; 

Row  2  will  correspond  to  Item  2  in  the  map;  etc. 

Columns  are  grouped  into  attribute  fields  and  each  attribute  field  must 
contain  a  single  type  of  data.  For  example.  Columns  1  throuc'ri  4   ay  be 
the  drill-hole  identifier,  columns  8  through  11  the  coal  thukne.  :3 ,  and  so 
forth.   The  length  of  an  attribute  field  must  be  the  same  within  each  row. 

Here  is  an  example  of  a  multiple  attribute  file  for  3  drill  hciles: 

149.0  5958 
lSi.5  5979 
220.5   5997 


DH-1 

2 

11.5 

136.5 

1.0 

DH-2 

6 

22.  0 

101.0 

73.6 

DH-3 

3 

36.8 

120.0 

80.2 

FIGURE  2.  LISTING  OF  THE  IMPORT  COMMAND. 


ENTER  COMMAND  ?   IMPORT  COALPTS .  IMP 
What  do  you  wish  to  call  the  new  map  ? 


:  COALPTS 
Enter  NAME  of  map  to  use  as  a  template  for  the  new  map  header  or  enter 
CARRIAGE  RETURN  to  start  map  header  from  scratch  : 

ENTER  SOURCE  OF  MAP   [  ] 

:  IMPORT 
ENTER   CREATION  DATE   [         ] 

:  1989 
ENTER   STUDYAREA  NAME  [         ] 

:  OZ 
ENTER   DESCRIPTION 
[  ] 

:  COAL  POINTS  FOR  TRAINING 
ENTER   MAP  VINTAGE   [     ] 

:  1989 
ENTER  NUMBER  OF  SUBJECTS  [     1] 

:  45 
ENTER  COORDINATE  SCALE  FACTOR  [     100.0000] 

:  100 

MAP  PROJECTION 

#  0  PROJECTION  IS  GEOGRAPHIC (LON/LAT) 
COORDINATE  UNITS  ARE:  DEGREES 

DO  YOU  WISH  TO  CHANGE  THE  PROJECTION  DESCRIPTION  [N] : 

Y 

PROJECTION (0-2  0)  ? 

1 
ELLIPSOID(0-19)  ? 

0 
LONGITUDE  OF  ANY  POINT  WITHIN  THE  UTM  ZONE  ? 

-108 
LATIT"DE  OF  ANY  POINT  WITHIN  UTM  ZONE    13  ? 

36 
IS  THIS  HEADER  INFORMATION  CORRECT  [Y] 

:  Y 
ENTER  DATAT-i^E 

1  =  POINT  11  =  (X,Y,Z)  POINT 

2  =  LINE  12  =  (X,Y,Z)  LINE 

3  =  POLYGON  13  =  (X,Y,Z)  POLYGON 
5  =  SAMPLE  ELEVATION  POINT 

:  1 

EXECUTING. . .PLEASE  WAIT 


IMPORT  COMPLETE  FOR  THE  NEW  MAP:  COALPTS 
45  ITEMS  AND     45  SUBJECTS  IN  THE  NEW  MAP 
THE  INPUT  FILE  IS  COALPTS . IMP 
DO  YOU  WISH  TO  DELETE  THE  INPUT  FILE   [N] 

:  N 

****  STOP 


FIGURE  3.   Listing  of  the  PROJECTION  command. 


ENTER  COMMAND  ?   PROJECTION  1 

What  do  you  wish  to  call  the  new  map  ?   PROJ.MAP 

WHAT  IS  THE  PROJECTION  OF  THE  INPUT  MAP  ? 

MAP  PROJECTION 
#  0  PROJECTION  IS  GEOGRAPHIC (LON/LAT)  COORDINATE  UNITS  ARE: 


DEGREES 


MAP  PROJECTION  MENU 

0  GEOGRAPHIC (LAT/LON) 

1  UNIVERSAL  TRANSVERSE  MERCATOR  11 

2  STATE  PLANE  12 

3  ALBERS  CONICAL  EQUAL  AREA  13 

4  LAMBERT  CONFORMAL  CONIC  14 

5  MERCATOR  15 

6  POLAR  STEREOGRAPHIC  16 

7  POLYCONIC  17 

8  EQUIDISTANT  CONIC  A  &  B  18 

9  TRANSVERSE  MERCATOR  19 

10  STEREOGRAPHIC  2  0 


LAMBERT  AZIMUTHAL  EQUAL  AREA 
AZIMUTHAL  EQUIDISTANT 
GNOMONIC 
ORTHOGRAPHIC 

VERTICAL  NEAR  SIDE  PERSPECTIVE 
SINUSOIDAL 

EQUIRECTANGULAR  (PLATE  CARREE) 
MILLER  CYLINDRICAL 
VAN  DER  GRINTEN  I 
OBLIQUE  MERCATOR  (HOTINE) 


PROJECTION (0-2  0)  ? 
0 

MAP  PROJECTION 
#  0  PROJECTION  IS  GEOGRAPHIC (LON/LAT) 

COORDINATE  UNITS  ARE:  DEGREES 
Press  CR  to  continue? 

WHAT  IS  THE  PROJECTION  FOR  THE  OUTPUT  MAP  ? 


MAP  PROJECTION 

0  GEOGRAPHIC (LAT/LON) 

1  UNIVERSAL  TRANSVERSE  MERCATOR   11 

2  STATE  PLANE  12 

3  ALBERS  CONICAL  EQUAL  AREA  13 

4  LAMBERT  CONFORMAL  CONIC  14 

5  MERCATOR  15 

6  rOLAR  STEREOGRAPHIC  16 

7  POLYCONIC  17 

8  EQUIDISTANT  CONIC  A  &  B  18 

9  TRANSVERSE  MERCATOR  19 

10  STEREOGRAPHIC  2  0 

PROJECTION (0-20)  ?   1 


MENU 

LAMBER1  AZIMUTHAL  EQUAL  AREA 
AZIMUTHAL  EQUIDISTANT 
GNOMONIC 
ORTHOGRAPHIC 

VERTICAL  NEAR  SIDE  PERSPECTIVE 
SINUSOIDAL 

EQUIRECTANGULAR  (PLATE  CARREE) 
MILLER  CYLINDRICAL 
VAN  DER  GRINTEN  I 
OBLIQUE  MERCATOR  (HOTINE) 


ELLIPSOID  MENU 


0 

CLARKE  1866  (default) 

10 

MODIFIED  EVEREST 

1 

CLARKE  1880 

11 

MODIFIED  AIRY 

2 

BESSEL 

12 

WALBECK 

3 

NEW  INTERNATIONAL  1967 

13 

SOUTHEAST  ASIA 

4 

INTERNATIONAL  1909 

14 

AUSTRAILIAN  NATIONAL 

5 

WGS  72 

15 

KRASOVSKY 

6 

EVEREST 

16 

HOUGH 

7 

WGS  66 

17 

MERCURY  19  60 

8 

GRS  198  0 

18 

MODIFIED  MERCURY  lb  08 

9 

AIRY 

19 

SPHERE  OF  RADIUS  637099"; 

ELLIPSOID (0-19)  ?    0 

LONGITUDE  OF  ANY  POINT  WITHIN  THE  UTM  ZONE  ?   -108 

LATITUDE  OF  ANY  POINT  WITHIN  UTM  ZONE    13  ?     3  6 


MAP  PROJECTION 


#  1  UNIVERSAL  TRANSVERSE  MERCATOR 

COORDINATE  UNITS  ARE:  METERS 
VALUES  OF  ENTERED  PARAMETERS 

ELLIPSOID  SELECTED  IS  TYPE 

SEMI -MAJOR  AXIS  OF  ELLIPSOID 

ECCENTRICITY  SQUARED  OF  ELLIPSOID 

ZONE 

LONGITUDE 

LATITUDE 


6378206 
0 

-108 
36 

YES 


4000000 
0067687 
13 
0000000 
0000000 


Output  projection  satisfactory  (Y/N)  [CR=Y] 

ENTER  SOURCE  OF  MAP:     IMPORT 

ENTER   CREATION  DATE:      1989 

ENTER   DESCRIPTION:      COAL  MAP  IN  UTM  PROJECTION 

ENTER  COORDINATE  SCALE  FACTOR:    100 


IS  THIS  HEADER  INFORMATION  CORRECT:   YES 


PROJECTING. . .PLEASE  WAIT 
****  STOP 


or  it  could  look  like  this  without  the  spaces  between  attribute  fields; 


DH-1211.5136.5  1.0149.05958 
DH-2 62 2. 010 1.073, 6 19 6. 65979 
DH-3 3 3 6. 8 12 0.080. 222 0.55997 


FACTS  TO  REMEMBER: 

There  may  be  up  to  2  00  attributes  per  drill  hole. 

Each  attribute  field  must  contain  a  single  type  of  data. 

Attribute  fields  for  each  drill  hole  must  be  of  the  same  length 

and  type. 

An  attribute  field  must  not  be  greater  than  7  0  characters 
in  length. 

The  order  of  the  rows  must  exactly  match  the  item  number  order  of  the 
drillhole-location  map.   (If  the  data  file  is  not  in  the  proper  order, 
use  the  UTILITY-ATTRIBUTE-RESEQUENCE  command  to  reorder  the  multiple 
attribute  file.) 


There  are  three  DATA  TYPES; 

CHARACTER  -  all  the  printable  ASCII  characters  including 
numerics,  symbols  and  spaces.    Characters  are  left 
justified  within  an  attribute  field. 
EXAMPLE:    DH-1 

REAL  -  positive  and  negative  real  numbers  as  large  as  1.8  x  10^^ 
and  as  small  as  1.0  x  10~63.   The  format  for  real  numbers  is 
deterimined  by  the  format  statement.   For  example,  F7 . 4  means 
there  are  a  total  of  7  spacss  with  4  numbers  to  the  right  of  the 
decimal . 

EXAMPLE:   44.9876 

INTEGER  -  positive  and  negative  integer  numbers  as  large  as  32,000. 
No  decimal  points  allowed.   An  integer  larger  than  3  2,000  must 
be  encoded  as  character  data  or  have  a  scaling  factor  a-.:  ied 
to  the  values.  An  inte  -r  is  right  justified  within  an 
attribute  field. 
EXAMPLE:   6678 


STEP  2.  CREATE  THE  DEFINITION  FILE 

The  user  must  create  a  definition  file  before  adding  the  attributes  to 
the  drillhole-location  map.  The  definition  file  stores  information  about 
each  attribute  field  within  a  data  file,  such  as: 

ATTRIBUTE  (KEY)  NAME   [up  to  10  characters  with  no  spaces] 
DESCRIPTION  of  the  attribute   [up  to  60  characters] 
DATA  TYPE   [integer,  real  or  character] 

FIELD  WIDTH  [number  of  spaces  for  the  attribute;  limitations 

according  to  data  type] 

*  A  DEFINITION  FILE  IS  UNIQUE  TO  THE  DATA  FILE  IT  IS  DEFINING.,  * 

Naming  convention  gives  the  suffix  .DEF  to  a  definition  file.   Here  are 
some  examples  of  definition  filenames: 

COAL. DEF        PI. DEF 

This  file  can  be  created  interactively  using  the  UTILITY-ATTDESCRIBE 
command.   However,  before  using  ATTDESCRIBE ,  determine  the  format  of  the 
attribute  file.   For  example,  the  first  line  in  a  file  COAL. DAT  might 
look  like  this: 

DH-1  2  11.5  136.5  1.0  149.0  5958 

Drill-hole  identification  (DH-1)  is  in  columns  1-4  and  is  character 

type. 
Number  of  coal  beds  (2)  is  in  spaces  5-6  and  is  an  integer. 
Total  coal  thickness  (11.5)  is  in  spaces  7-11  and  is  a  real  number 

with  one  place  to  the  right  of  the  decimal. 
Overburden  (136.5)  is  in  spaces  12-17  and  is  a  real  number  with  one 

place  to  the  right  of  the  decimal. 
Total  interburden  (1.0)  is  a  real  number,  spaces  18-21. 
Total  depth  (149.0)  is  a  real  number,  spaces  22-27. 
Surface  elevation  (5958)  is  an  integer,  spaces  28-32. 

Note  that  every  space  has  been  accounted  for,  as  ATTDESCRIBE  assume-^ 
there  are  no  spaces  between  attribute  fields.   If  the  file  has  spa  "ss 
between  fields,  either  edit  these  spaces  out  of  the  attribute  file  or 
create  the  definition  file  and  then  edit  the  definition  file. 

Refer  to   £GURE  4  for  an  example  of  the  ATTDESCRIBE  command. 

STEP  3  ■  -  -DP  THE  ATTRIBUTES 

FIGURE  5  is  an  example  of  the  UTILITY  -  ADD  ATTRIBUTE  command,  which 
adds  the  attribures  to  the  map  by  using  the  definition  file  (.DEF)  and 
the  data  file  (.DAT  or  .ATT).    This  command  creates  a  binary  file  that 
has  a  mapname.AT  filename.   For  example,  a  map  named  COAL.UTM  would  have 
a  binary  attribute  file  named  COAL.UTM. AT. 


FIGURE   4.      Listing  of  the  UTIUTY-ATTDESCRIBE  conunand. 


ENTER 
ENTER 

1  = 

2  = 

-3     — 
^     — 

4   = 

b  = 

7  = 

8  = 

9  = 

10  = 

11  = 

12  = 

13  = 

14  = 

15  = 

16  = 


COMMAND  ?   UTIL 
MOSS  UTILITY  OPTION 
TERMINATE  UTILITY  SESSION 


[DEFAULT] 


DATABTEST 

ATTRIBUTE 

ATTDES 

SUB2AT 

APROJ 

BROWZ 

PLOT . LEGEND 

MAKE . LOGO 

SUBEDIT 

SET. LEVEL 

TRANSFORM 

QUAD 

DLG3 

MAPIDX 

XYSUBJECT 


(MOSS  MAP  NAMES  SUPPORT) 

(MOSS  MULTIPLE  ATTRIBUTE  SUPPORT) 

(BUILD  MULTIPLE  ATTRIBUTE  DEFINITION  FILE) 

(SUBJECT  TO  MULTIPLE  ATTRIBUTE  INPUT) 

(MOSS  MAP  NAMES  PROJECTION  ASSIGNMENT) 

(MOSS  MAP  NAMES  HEADER  LISTING) 

(BUILD  PLOTTER  LEGEND  FILE) 

(BUILD  PLOTTER  LOGO  FILE) 

(MAP  SUBJECT  EDIT  PROGRAM) 

(BUILD  POLYCELL  TRANSLATION  FILE) 

(TRANSFORM  COORDINATES  TO  A  PROJECTION) 

(MAKE  A  QUAD  MAP  IN  IMPORT/EXPORT  FORMAT) 

(USGS  DLG  ASCII  TO  MOSS) 

(MAKE  INDEX  MAP  OF  PROJECT  IN  IMPORT/EXPORT  FORMAT) 

(REFORMAT  POINT  DATA  TO  MOSS) 


***  ATTRIBUTE  DESCRIPTION  PROGRAM  *** 

What  do  you  wish  to  call  the  Definition  File? 

COAL.DEF 
Provide  a  ten  character  description  for  "KEY" 

DRILLHOLE 
Provide  a  6  0  character  description  of  this  "KEY" 

DRILL  HOLE  NUMBER 
The  field  type  of  this  "KEY"  is 
1  -  Integer   2  -  Floating  point   3  -  Character 
Select:   3 
What  is  the  field  length  of  this  "KEY"?    6 

Attribute:         1    Key:       DRILLHOLE 

Description:   DRILL  HOLE  NUMBER 

Type:      CHARACTER     Length:        6 


1  [CR=Exit] 


Edit  options; 


[0]  -  Proceed  to  next  attribute  (# 

1  -  Edit  attribute  KEY 

2  -  Edit  attribute  DESCRIPTION 

3  -  Edit  field  TYPE  and  LENGTH 

4  -  QUIT 

5  -  ABORT 


2) 


Select:  0 


Provide  a  ten  character  description  for  "KEY" 

THICKNESS 
Provide  a  60  character  description  of  this  "KEY" 

THICKNESS  OF  THE  COAL  BED 
The  field  type  of  this  "KEY"  is 


2  [CR=Exit] 


1  -  Integer  2  -  Floating  point  3  -  Character 

Select!   2 

What  is  the  field  length  of  this  "KEY"?   4 

How  many  digits  to  the  right  of  the  decimal  point?  1 


Attribute: 

Description: 

Type: 

Edit  options; 


2     Key:       THICKNESS 
THICKNESS  OF  THE  COAL  BED 
REAL    Length:       4.  1 

[0]  -  Proceed  to  next  attribute  (# 

1  -  Edit  attribute  KEY 

2  -  Edit  attribute  DESCRIPTION 

3  -  Edit  field  TYPE  and  LENGTH 

4  -  QUIT 

5  -  ABORT 


3) 


Select: 


Provide  a  ten  character  description  for  "KEY"    3  [CR=Exit] 

STRUC.TOP 
Provide  a  60  character  description  of  this  "KEY" 

ELEVATION  OF  TOP  OF  COAL  BED 
The  field  type  of  this  "KEY"  is 
1  -  Integer   2  -  Floating  point   3  -  Character 
Select :   1 
;vhat  is  the  field  length  of  this  "KEY"?   5 

Attribute:  3  Key:  STRUC.TOP 
Description:  ELEVATION  OF  TOP  OF  COAL  BED 
Type:        INTEGER    Length:        5 


Edit  options; 


[0]  -  Proceed  to  next  attribute  (# 

1  -  Edit  attribute  KEY 

2  -  Edit  attribute  DESCRIPTION 

3  -  Edit  field  TYPE  and  LENGTH 

4  -  QUIT 

5  -  ABORT 


4) 


Select:  4 


****  STOP 


1 

= 

TERMINATE  U 

2 

= 

DATABTEST 

3 

= 

ATTRIBUTE 

4 

= 

ATTDES 

5 

= 

SUB2AT 

6 

= 

APROJ 

7 

= 

BROWZ 

8 

= 

PLOT. LEGEND 

9 

= 

MAKE . LOGO 

10 

= 

SUBEDIT 

11 

= 

SET. LEVEL 

12 

= 

TRANSFORM 

13 

= 

QUAD 

14 

= 

DLG3 

15 

= 

MAPIDX 

16 

= 

XYSUBJECT 

FIGURE  5.   Sample  run  of  the  MOSS  command  unUTY  to 

add  attributes  to  the  drillhole-location  map. 

ENTER  COMMAND  ?   UTILITY 

ENTER  MOSS  UTILITY  OPTION 

i-ILITY  SESSION    [DEFAULT] 
(MOSS  MAP  NAMES  SUPPORT) 
(MOSS  MULTIPLE  ATTRIBUTE  SUPPORT) 
(BUILD  MULTIPLE  ATTRIBUTE  DEFINITION  FILE) 
(SUBJECT  TO  MULTIPLE  ATTRIBUTE  INPUT) 
(MOSS  MAP  NAMES  PROJECTION  ASSIGNMENT) 
(MOSS  MAP  NAMES  HEADER  LISTING) 
(BUILD  PLOTTER  LEGEND  FILE) 
(BUILD  PLOTTER  LOGO  FILE) 
(MAP  SUBJECT  EDIT  PROGRAM) 
(BUILD  POLYCELL  TRANSLATION  FILE) 
(TRANSFORM  COORDINATES  TO  A  PROJECTION) 
(MAKE  A  QUAD  MAP  IN  IMPORT/EXPORT  FORMAT) 
(USGS  DLG  ASCII  TO  MOSS) 

(MAKE  INDEX  MAP  OF  PROJECT  IN  IMPORT/EXPORT  FORMAT) 
(REFORMAT  POINT  DATA  TO  MOSS) 
:  3 

PLEASE  ENTER  MAP  NAME:    COAL.PTS 

THE  NUMBER  OF  ATTRIBUTES  IS      0 
THE  NUMBER  OF  ITEMS  IS  45 

PLEASE  ENTER  DESIRED  OPTION 

1  =  EXIT    [DEFAULT] 

2  =  ADD  A  NEW  ATTRIBUTE 

3  =  UPDATE  AN  EXISTING  ATTRIBUTE 

4  =  CHANGE/DELETE  KEY  OR  DESCRIPTOR 

5  =  LIST  ATTRIBUTE  FIELDS 

6  =  SEARCH  AN  ATTRIBUTE  FIELD 

7  =  RESEQUENCE  INPUT  DATA  FILE 

8  =  DELETE  THE  ATTRIBUTE  FILE 
:  2 

will  you  be  using  a  definition  file?  [cr  =  no]:     yes 
enter  the  name  of  the  attribute  defintion  file:     coal.def 
please  enter  the  name  of  the  input  data  file:      coal. dat 

the  number  of  attributes  is    5 
the  number  of  items  is        45 

plea::"  enter  desired  option 

1  exit   [default] 

2  add  a  new  attribute 

3  -  update  an  existing  attribute 

4  =  change/delete  key  or  descriptor 

5  ==  list  attribute  fields 

6  =  search  an  attribute  field 

7  =  resequence  input  data  file 

8  =  delete  the  attribute  file 

:  1 
****  STOP 
MOSS  UTILITY  SESSION  COMPLETED 


Now  the  user  can  use  the  following  commands: 

DESCRIBE  mapname  ATTRIBUTE  to  view  the  attribute  file. 
SELECT,  WINDOW  and  PLOT  the  drillhole-location  map. 
LEGEND  mapid  LABEL  to  plot  the  subject  or  an  attribute 

next  to  the  drill-hole  locations. 
SELECT  mapname  ATTRIBUTE  to  select  a  specific  attribute  for 

plotting  on  the  screen. 
REPORT  to  create  a  table  that  lists  the  attributes  as 

specified  by  the  user. 
QUERY  ATTRIBUTE  to  query  a  specific  point  on  the  plotted 

drillhole-location  map. 
BSEARCH  to  perform  complex  Boolean  retrievals  of  information  from 

the  attribute  file  of  a  map. 

If  necessary,  attributes  can  be  edited  with  the  UTILITY-ATTRIBUTES 
command  (Option  3-update  an  existing  attribute)  or  the  EDITATT  command. 
The  first  command  will  use  a  data  file  to  modify  all  the  values  for  an 
existing  attribute;  the  second  command  uses  the  crosshairs  to  locate  a 
specific  drill  hole  and  modify,  by  key  entry,  a  specific  attribute,  or  all 
the  attributes,  for  that  drill  hole. 

SECTION  FIVE:  CREATE  A  BORDER  MAP  FOR  THE  STUDY  AREA 

Before  generating  isopach  and  structure  maps,  a  border  map  must  be 
created  that  will  be  used  as  the  "window"  for  all  MAPS  analyses.   The 
border  map  is  essential  in  the  analyses  as  all  the  cell  maps  must  have 
the  same  number  of  cells  (rows  and  columns)  and  the  same  cell  size. 

There  are  several  ways  to  do  this  depending  on  what  maps  are  available. 

If  there  is  a  specific  lease/tract  boundary  for  the  study  area: 

SELECT,  WINDOW  and  PLOT  the  landlines  (township,  range)  map(s). 

WINDOW  the  active  ID(s)  of  the  landline  map(s). 

SELECT  and  PLOT  the  lease  boundary  that  is  within  the  landline 

boundary. 
Use  the  GENERATE  command  and  the  crosshairs  to  create  a  rectangle 

just  oustide  the  lease  boundary.   Name  the  new  map  BORDER. 
SELECT  and  WINDOW  the  BORDER  map.   Use  this  window  for  all 

MAPS  analyses  for  a  project. 

If  there  is  no  specific  lease  or  other  boundary,  use  the  drillhci--- 
location  map  to  creare  a  BORDER  map: 

SELECT,  WINDOW  and  PLOT  the  location  map. 
Use  the  command  ZOOM  POINT  with  a  factor  of  .9  or  .3. 
ERASE  the  screen  and  rePLOT  the  location  map. 

Use  the  GENERATE  command  to  create  a  rectangle  by  positioning  the 
crosshairs  just  inside  the  corners  of  the  enlarged  map.   Name 


the  new  map  BORDER. 
SELECT  and  WINDOW  the  BORDER  map.  Use  this  window  for  all  the 
MAPS  analyses  for  the  project. 

SECTION  SIX:  USE  TESTGRID  TO  DETERMINE  THE  CELL  SIZE  FOR  THE  GRID 

Prior  to  the  gridding  process,  the  user  must  determine  an  appropriate 
grid  cell  size.   The  TESTGRID  command  can  be  used  to  visually  determine 
what  cell  size  would  be  most  useful. 

PLOT  the  drillhole-location  map 
Use  the  TESTGRID  command 

A  general  "rule" —  no  more  than  1  point  to  a  cell;  points  not 
more  than  10  cells  away  from  each  other. 

SECTION  SEVEN:  THE  ABC'S  OF  GRIDDING 

The  MOSS  GRID  command  converts  a  point-value  map  (e.g.,  observed 
thickness  of  coal)  into  a  uniform  set  of  rows  and  columns,  thus  defining 
a  grid  made  up  of  "cells".   Each  cell  is  assigned  a  specific  value  that 
has  been  interpolated  by  MOSS  from  the  observed  point -values.   The  grid 
can  then  be  used  to  generate  a  contour  map  (e.g.,  isopach  map)  or  in 
further  analyses  (.e.g.,  adding  grids  of  several  coal  beds  to  obtain  a 
grid  of  total  coal  thickness) . 

GRIDDING  ALGORITHMS  IN  MOSS; 

There  are  5  gridding  algorithms  available  in  MOSS.   Each  one  treats 
the  gridding  process  differently,  and  the  user  should  select  an  algorithm 
depending  on  the  configuration  of  the  data  (e.g.,  clustered,  widely 
scattered,  etc.)  and  his/her  geostatistical  knowledge.   The  gridding 
algorithms  available  in  MOSS  are  as  follows: 

4-POINT/QUADRANT  WEIGHTED  AVERAGE:   Searches  in  four  directions 
(quadrant)  from  a  target  cell  to  find  the  four  nearest  observed  values 
and  then  calculates  the  weighted  average  of  those  values  for  the  cell. 
The  value  (z)  for  a  cell  is  determined  by  this  equation: 

z   =   zi  /di   +   22  7^2   +   Z3  7^3   +  24  /d4 

1/di  +   l/d2   +  l/d3   +  1  'd4 

where: 

z      =  the  cell  value  to  be  determined 
1-..4  =  the  value  of  the  observed  points 
d     =  the  distance  from  an  observed  point  to  the 
center  of  the  target  cell 


This  algorithm  works  better  than  others  when  there  is  a  clustering  of  the 

data. 

8-POINTS  WEIGHTED  AVERAGE :   Searches  in  all  directions  for  the  8  nearest 

observed  point -values  from  the  target  cell-   The  z  value  of  the  cell  is 
then  calculated  by  a  similar  formula  as  above  except  for  8  points  rather 
than  4 .  The  8-POINT  algorithm  produces  a  smoother  contour  map  than  the  4- 
POINT. 

KRIGiNG:  Estimates  the  value  of  a  spatially  distributed  variable  from 
adjacent  values  while  considering  their  inter independence  as  shown  in  a 
variogram.   The  user  should  not  attempt  to  run  the  KRIGING  option  without 
some  understanding  of  the  theoi-y  of  geostatistics. 

QUINTIC  SPLINE  INTERPOLATION:   Divides  the  area  into  triangles  using  the 
points  as  vertices  of  the  triangles.   It  then  defines  a  5th  degree 
(quintic)  polynomial  over  each  triangle  in  such  a  manner  that  the 
polynomials  of  adjacent  triangles  match  smoothly  (spline) .   This  gridding 
algorithm  is  good  for  interpolating  within  the  convex  hull  of  the  data 
points;  not  good  for  outside  this  hull. 

NOTE:  The  QUINTIC  SPLINE  method  cannot  be  used  if  two  points 
have  the  same  geographic  coordinates.   The  VARIOGRAM  command  can  be 
used  to  identify  such  points  as  it  lists  the  item  niimbers  of  the 
duplicate  points.  The  SELECT  command  can  then  be  used  to  not-select 
the  duplicate  item  numbers. 

NEAREST  NEIGHBOR  (POLYGON):    Searches  for  the  value  of  the  closest  data 
point  to  the  target  cell  and  assigns  this  value  to  the  cell.   If  more 
than  one  data  point  is  available  to  the  target  cell,  it  is  assigned  the 
average  of  the  data  points.   The  resulting  cell  map  appears  to  made  of 

polygons. 

Several  of  the  above  gridding  algorithms  will  ask  for  the  size  of  a 
"roving  window"  (a  moving  matrix  of  cells) .   The  user-specified  window 
size  defines  the  number  of  cells  (in  both  x-  and  y-dimensions)  that  the 
gridding  algorithm  will  search  for  the  required  number  of  observed  point- 
values.   If  three  points  cannot  be  found  within  the  window,  the  target 
cell  is  given  the  value  of  zero. 

SECTION  EIGHT:  CREATE  AN  ISOPACH-GRiD  MAP 

STEP  1.   SELECT  the  drillhole-location  map  by  the  thickness  att,:  bute  or 
BESEARCH  the  map  using  the  Boolean  operators,  such  as  in  the  fox_owing 
example: 

BSEARCH  COAL.PTS 

THICKNESS  GT  0 


STEP  2 .   Use  the  GRID  command: 

WINDOW  the  BORDER  map   [as  created  in  SECTION  FIVE] 
GRID  the  thicJoiess  map   [using  cell  size  determined  in 

SECTION  SIX] 

An  example  of  the  QUINTIC  SPUNE  gridding  method  is  as  follows: 

GRID  the  active  table  ID  #  for  thiclcness  values 
MENU  #?  5  [Quintic  Spline] 

MASKING  MAP?  NO 
SCALER:  1 

ACRES  OR  METERS?  M 
CELL  WIDTH?  500 
CELL  HEIGHT?  500 
REENTER  PARAMETERS?   N 

NO.  OF  CLOSEST  DATA  POINT  TO  USE?    7 
IS  THIS  A  RESTART?    N 

The  user  should  create  and  plot  various  types  of  grids  and  choose  the 
best  one  for  the  data  in  the  study  area.   The  QUINTIC  SPUNE  algorithm 
creates  the  grid  faster  but  needs  a  wide-spread  control  by  the  data  or  it 
will  see  trends  where  the  data  are  absent.   For  clustered  data,  use  the 
4-POINT/QUADRANT  WEIGHTED  algorithm. 

The  grid  may  then  be  SELECTed  and  PLOTted  on  the  screen.   The  drillhole- 
location  map  and  other  available  maps  (e.g.,  lease  boundary)  can  be 
PLOTted  over  the  isopach-grid  map.  See  FIGURE  6  CA  through  C)  for 
examples  of  the  4-POINT,  8-POINT  and  QUINTIC  SPUNE  grids  using  the  same 
data  set. 

SECTION  NINE:  CREATE  A  STRUCTURE-GRID  MAP 

The  process  is  the  same  as  the  isopach  map  except  that  the  attribute  for 
the  elevation  for  the  top  of  the  coal  bed  is  used  instead  of  the  coal 
thickness.   See  FIGURE  7  for  an  example  of  a  structure  map  generated 
by  the  QUINTIC  SPUNE  method. 

Other  MOSS  commands  that  can  be  used  with  the  grid  map: 

-  The  SCAN  command  in  MAPS  can  be  used  on  the  grid-ma:  tc  el:..aj  -v-r.e 
excessive  contour  islands  in  the  plot. 

-  The  CONTOUR  command  in  MOSS  will  create  a  line  map  from  tr-e  c  ;.d 
map.   The  SMOOTH  command  in  MOSS  will  make  the  contour  line  map 
smoother  each  successive  time  it  is  used. 

-  The  contour  lines  can  be  labeled  using  the  LEGEND  -  CONTOUR 
command  or  the  AUTOLABEL  command. 

-  Other  maps  can  be  PLOTted  over  the  contour  map,  such  as  the 
drillhole-location  map,  lease  map,  landlines,  roads,  streams,  etc. 


FIGURE   6.    Isopach  map  generated  by  the    following   gridding   algorithms: 
^     A-POIKT"  OllADRANT:     (B)    a-POINT  WEIGHTED ?    AND    (C)     QUINTIC  SPUNE. 


FIGURE   7 .      Structure  map   generated  by  the  QUINTIC  SPLINE  gridding 


-  For  a  hardcopy  map,  the  maps  can  be  plotted  to  a  specified  scale 
using  the  PENPLOT  command. 

NOTE:   The  SCAN  and  SMOOTH  commands  will  often  smooth  out  contours  so 
that  they  no  longer  honor  the  data  points,  i.e.,  the  10'  contour  may  not 
go  through  a  10'  data  point.   Remember,  however,  that  it  is  the  grid  (not 
the  contour  line  map)  used  in  MAPS  analyses,  such  as  calculating  coal 
tonnage. 

SECTION  TEN:  EDITING  THE  GRID.  IF  NECESSARY 

Once  the  grid  is  generated,  the  geologist/engineer  can  determine  if  it 
geologically  acceptable.   If  not,  either  try  other  parameters  (cell  size, 
window)  or  another  gridding  algorithm. . .or  edit  the  grid. 

Editing  may  consist  of  one  or  more  of  the  following  methods: 

1.  Use  the  DESCRIBE  command  to  see  if  the  grid  has  cell  values  that 
are  negative  or  are  higher  than  the  the  largest  observed  data  value. 
This  would  be  especially  true  if  the  QUINTIC  SPLINE  algorithm  was 
used  as  it  creates  trends  where  observed  data  is  sparse  or  lacking. 
The  RENUMBER  command  can  be  used  to  change  these  values  in  the  grid. 
Here  is  an  example: 

MAPS 

DESCRIBE  THICK.  GRD  HEADER 
WINDOW  BORDER 

RENUMBER  THICK. GRD  FOR  NEWTHK.GRD  ASSIGN  0  TO  -50  TH  0  , 
ASSIGN  25  TO  25  TH  50 

2.  Data  errors  can  sometimes  be  spotted  when  plotting  the  grid  and 
there  are  isolated  areas  with  tight  concentric  contours.   These 
data  mirrors  can  also  be  seen  as  spikes  (positive  or  negative)  in  a 
three-dimensional  display  of  the  grid  using  the  3D  command  in  MAPS, 
The  data  errors  can  be  deleted  by  one  of  2  methods: 

a.  PLOT  and  QUERY  the  map  to  determine  the  item  numbers  of  the 
drill  holes  that  have  data  errors.   Then  SELECT  the 
drillhole-location  map  by  ITEM  and  use  the  non-select  option 

'a  backslash) . 

SELECT  mapname  ITEM 

\45  57  [not  items  45  and  57] 

This  selection  can  then  be  saved  for  a  new  map  using  the  SAVE 
command.  The  gridding  process  must  then  be  redone  using  the 
new  map. 

b.  Use  the  EDMAP  command  to  delete  the  points  using  the 


crosshairs  and  create  a  newmap.  Again  the  gridding  must  be 
redone  using  the  new  map. 

3.  Use  the  EDMAP  command  to  add  control  points.  Use  minimum 
number  of  points  and  reGRID  before  adding  more  points.   Control 
points  work  well  in  the  QUINTIC  SPUNE  method,  but  often  just  add 
more  circles  to  the  8-POINT  grid. 

4.  Change  individual  or  blocks  of  cell  values  using  the  CUT 
command. 

5.  CONTOUR  the  grid  and  then  use  EDMAP  to  change  the  contour  lines. 
However,  this  only  changes  the  contour  map  and  not  the  gridi   To 
get  the  map  back  into  grid  format,  use  the  CTOG  command  and. create 
a  grid  from  the  corrected  contour  map. 

When  the  grid  is  geologically  acceptable  to  the  geologist/engineer,  it 
can  be  used  to  generate  plots  (isopach  maps,  structure  maps)  or  in 
further  data  analyses  such  as: 

MATH  DEMgrid  -  structure . grid  FOR  overburden . grid 

MATH  isopachl.grid  +  isopach2 .grid  FOR  total. coal. grid 

SECTION  ELEVEN:  COAL  RESERVE  CALCULATIONS 

In  calculating  coal  reserves,  the  thickness  values  must  first  be  limited 
to  the  lease  (or  sections  within  a  lease)  boundary  so  that  only  the  coal 
within  that  boundary  is  considered  in  the  calculations.   The  boundary, 
which  is  a  vector  map,  must  be  rasterized  to  a  create  a  cell  map  that 
has  the  same  cell  size  and  rows  and  columns  as  the  thickness-grid  map. 
The  thickness-grid  (feet  of  coal  within  the  lease  boundary)  and  the 
rasterized  boundary  map  (acres)  are  used  in  the  tonnage  calculations  to 
obtain  acre-feet.   A  coal  factor  can  also  be  used.   The  steps  to  do  this 
are  as  follows: 

STEP    1.  USE  THE  LPOVER  COMMAND  TO  LIMIT  CONTOURS  TO  LEASE  BOUNDARY 

SELECT  and  WINDOW  the  BORDER  map. 

SELECT  the  thickness  grid  and  CONTOUR  it  to  create  a  line  map 

(Type  2) . 
SELEC"  the  contoured  map  (line  map) . 
SELECT  the  lease  boundary  map. 
LPOVER  the  lease  and  contour  maps  to  create  a  new  map. 


STEP    2.       CREATE   A   CELL  MAP   OF   THE    LEASE/TI^CT   MAP 

FACTS  TO  REMEMBER: 

Use  the  BORDER  map  as  the  window  in  all  cases  so  that  the  resultant 
cell  map  has  the  same  number  of  rows  and  columns  as  the  thickness 

grid. 

Use  the  same  cell  size  for  all  grids  and  rasterized  maps  for  the 
the  same  study  area. 

If  the  tonnage  is  to  be  calculated  by  the  lease  boundary: 

Rasterize  the  boundary  map,  such  as  the  lease/tract  boundary  ■ 
map,  using  the  MOSS  command  POLYCELL  or  the  MAPS  command  RASTERIZE. 

SELECT  lease. map  all. 
WINDOW  the  BORDER  map. 
POLYCELL  the  lease. map. 

If  tonnage  by  section  within  the  lease  boundary  is  required,  then  the 
following  extra  steps  must  be  taken: 

SELECT  the  lease. map  and  the  landlines  map. 
WINDOW  the  BORDER  map. 

INTERSECT  the  landlines  map  with  the  lease  map. 
POLYCELL  the  intersected  map. 

STEP  3.   SET  LINES  PER  PAGE  FOR  THE  REPORT 

Before  running  the  TOTAL  command  (for  coal  tonnage)  ,  set  the  lines  per 
page  for  printer  output  to  100  by  typing  PAGE  100.   This  will  change  the 
number  of  lines  per  page  from  3  0  to  100  so  that  a  title  line  does  not 
appear  in  the  middle  of  the  a  tonnage  report - 

STEP  4 .  USE  TOTAL  COMMAND  TO  CALCULATE  TONNAGE 

Here  is  an  example  of  using  the  TOTAL  command  in  MAPS  to  calculate  coal 

tonnage. 

MAPS 

WINDOW  BORDER 

TOTAL  ^HICK.GRD  FACTOR  1770  BY  LEASE.  POL  LABEL  TONS  FOR  REI    i:. 

w\  re  THICK. GRD  in  the  thickness  grid  (or  the  NEWTHK.GRD  if 
renumbered) 
FACTOR  is  the  coal  factor,  in  this  case  1770 
LEASE. POL  is  the  rasterized  lease  map 
LABEL  TONS  replaces  the  default  label  (PRODUCT)  in 

the  report 
FOR  REPT.LE  is  the  name  given  the  output  report  file 
that  will  list  the  coal  tonnage  for  the  lease 


Or  if  tonnage  by  section  within  the  lease  boundary  is  required: 

MAPS 

WINDOW  BORDER 

TOTAL  THICK. GRD  FACTOR  1770  BY  SECTION. POL  LABEL  TONS  FOR  REPT.SEC 

STEP  5.   PRINT  THE  REPORT 

The  output  report  will  be  a  text  file  in  your  directory.   Exit  MAPS  and 
MOSS  and  print  the  generated  tonnage  reports  at  the  PRIMOS  level.   See 
FIGURE  9  for  an  example  of  a  coal  tonnage  report. 

NOTE:   Tonnage  within  a  specified  boundary  will  not  be  accurate  because 
the  acreage  changes  when  the  lea^i-^  map  is  rasterized.    This  error  is 
inherent  in  working  with  cell  maps.   The  user  can  obtain  the  correct 
acreage  for  the  lease  by  using  the  AREA  command  in  MOSS  and  then  hand 
calculate  the  tonnage  by  using  the  average  thickness  as  generated  by  the 
TOTAL  command  in  MAPS. 

SECTION  TWELVE:  USING  THE  DIGITAL  ELEVATION  MODEL  (DEM^  TO  CALCULATE  THE 
OVERBURDEN  AND  LOCATE  THE  PROBABLE  OUTCROP  LINE 

STEP  1.   LOADING  THE  OEM'S 

Call  your  State  GIS  Coordinator  to  obtain  both  the  DEM  tape  and  the 
position  numbers  for  the  7.5  minute  guadrangle(s)  required  for  the  study 
area.   (It  may  be  necessary  to  know  the  USGS  ID  for  the  quadrangles, 
e.g.,  B344106.)   The  system  operator  will  then  load  the  specified  tape  on 
the  PRIME  and  inform  the  user  of  the  drive  number. 

Then  type  in  the  following: 

LOGIN  to  the  PRIME 

DOWN  directory  [to  user  directory] 

ASSIGN  MT#  [#  =  no.  of  the  mount  drive  on  which 

MAGNET  the  tape  is  loaded] 

READ 

MTU  NO  .?  =  #  [#  =  no.  of  mount  drive] 

FILE  NO,?  =  #  [#  =  file  number  for  the  quadrangle] 

LOGICAL  RECORD  LENGTH?  =  1024 

BLOCKING  FACTOR?  =  4 

TYPE?   ASCII 

DISCFILE?  dam. quad         [Enter  name  to  call  the  output  DEM] 

Repeat  for  as  many  quadrangle  DEM  files  needed  for  the  study  area. 
When  a  '>'  prompt  appears,  type  in  QUIT  to  exit  MAGNET. 

UNASSIGN  MT#  -UNLOAD        [rewinds  tape;  #  =  no.  of  mount  drive ^ 
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FIGURE  9.      Listing  of  coal  tonnage  within  each  section  of  a  lease 


STEP  2.   CONVERTING  THE  DEM 

This  section  will  take  the  user  through  the  necessary  steps  to  prepare  a 
DEM  (Digital  Elevation  Model)  so  that  it  has  the  same  projection,  cell 
size,  map  size  (rows  and  columns),  and  unit  of  measurement  (feet  rather 
than  meters)  as  the  thickness  and  structure  grids  so  that  it  can  be  used 
in  MAPS  analyses.    For  example,  the  structure  grid  can  be  stibtracted 
(using  MATH  command)  from  the  DEM  to  obtain  the  overburden  grid. 

A.  Log  into  MOSS  and  then  into  MAPS. 

B.  IMPORT  the  DEM  file: 

IMPORT  dem.quad  FORMAT  DEM  7.5  TYPE  8  for  dem.map 

where  dem.quad  is  the  output  DEM  file  from  Step  One  and 

dem.map  is  name  for  the  output  map  generated  by  the 
IMPORT  command. 

C.  MERGE  the  DEM  files  if  more  than  one  was  needed  for  the  study  area: 

MERGE  dem.mapl  with  dem.map2  for  demproject 

D.  RePROJECT  the  DEM  file  into  projection  used  with  other  maps  in  the  study 
area: 

PROJECT  demproject  FOR  dem.rep 

E.  To  convert  to  the  cell  size  of  the  DEM  to  the  cell  size  of  the 
thickness  and  structure  grids,  use  the  RESAMPLE  command: 

RESAMPLE  dem.rep  FOR  dem.resam  SIZE  # 
where  #  is  the  cell  siz  ; 

F-  DESCRIBE  the  resampled  DEM  and  the  thickness  or  structure  grid  to 
see  if  the  number  of  rows  and  columns  are  the  same  in  each  map.   Most 
likely  the  DEM  will  be  a  larger  map  and  the  CUT  command  must  be  used  to 
reduce  the  DEM  grid  to  the  same  size  as  the  thickness  and  structure  grids 
in  order  to  do  any  analyses: 

MAPS 

WINDOW    the    BORDER   map. 

CUT  dem.resam  WINDOW  FOR  dem.cut 

G.  To  ccnv-rt  elevations  frcm  meters  to  feet  in  the  DEM  file: 

MATH  dem.cut  *  3.230839894  FOR  dem.ft 


STEP  3.  COMPARE  THE  DEW  ^fAP  TO  DRILL-HOLE  COLLAR  ELEVATIONS 

The  collar  elevation  at  the  drill-hole  location  should  be  checked  against 
the  surface  elevation  as  given  by  the  DEM  at  that  location.   If  there  is 
a  sxibstantial  discrepancy  between  the  collar  elevation  and  the  surface 
elevation  of  the  DEM,  these  differences  will  have  to  be  evaluaited  by  the 
geologist/ engineer  and  corrections  made,  if  necessary,  to  the  structure 
and  thickness  values.    The  procedure  to  do  this  is  as  follows: 

Enter  MAPS. 

WINDOW  the  BORDER  map. 

RASTERIZE  the  drillhole-location  map  by  the  collar  elevation 
attribute  for  a  TYPE  7  map  using  the  same  cell  height 
and  width  as  used  with  the  thickness  grid: 

RASTERIZE  DH.MAP  BY  ATTRIBUTE  #  FOR  COLLAR. RAS  , 
HEIGHT  500  WIDTH  500 

where , 

DH.MAP  is  the  drillhole-location  mapname 

#  is  the  id  number  of  the  collar  elevation  attribute 

COLLAR.  RAS  is  the  name  given  the  output  map 

500  is  the  same  cell  size  in  the  thickness  grid 

Subtract  the  rasterized  map  from  the  DEM: 

MATH  DEM. FT  -  COLLAR. RAS  for  DIFF 

TOTAL  the  resulting  map  by  the  rasterized  map: 

TOTAL  DIFF  BY  COLLAR. RAS  FOR  REPT.DIFF 

The  report  will  be  generated  in  the  user's  directory.   The  report  can 
then  be  listed  on  the  screen  or  printed  to  hard  copy.   The  column  in  the 
report  that  lists  the  total  is  the  difference  between  the  DEM  elevation 
at  a  specific  point  and  the  collar  elevation. 

SECTION  THIRTEEN:  GENERATE  AN  OVERBURDEN  MAP 

SUBTRACT  the  structure  grid  from  the  DEM 

Remember,  the  DEM  must  have  elevation  coverted  to  feet,  cell  size  ~he 
same  as  the  structure  grid  and  cut  to  the  same  window. 

See  FIGURE  10  for  an  example  of  a  generated  overburden  map. 
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SECTION  FOURTEEN:  GENERATE  A  MINE-STRIPPING  RATIO  MAP 

A  mine-stripping  ratio  map  can  be  generated  to  limit  the  extent  of  the 
overburden,  interburden  or  coal: 

POLYCELL  the  bed  limit  polygon 
RENUMBER  the  polygon  area  to  1 
RENUMBER  the  rest  of  the  map  to  0 

Using  the  MATH  command  add  the  overburden  and  interburden  grids  together 
and  add  the  coal  bed  grids  together  as  follows: 

MATH  (  overburden  +  interburden  )  /  (  coal^  +  coal 2  +  coal^  )  , 
*   factor  FOR  strip. map 

The  factor  is  determined  by  the  user.   An  example  of  the  factor 
is: 

1613.33  cubic  yds.  /  acre-foot =     1,03 

(1742  tons/acre- foot)  *  .90  recovery  factor 

NOTE:  If  overburden  includes  coal  beds,  then  use  MATH  to  subtract  these 
beds  from  the  total  overburden: 

MATH  (  overburden  -  sum  coal  thickness  +  interburden)  /  , 
sum  of  coal  thickness  *  factor  FOR  strip. map 

See  FIGURE  11  for  an  example  of  a  generated  mine-stripping  ratio  map. 

SECTION  FIFTEEN:  CREATING  ISOPLETH  MAPS 

Isopleth  maps  can  be  generated  by  SELECTing  the  point  map  by  a  chemical 
attribute,  such  as  ash,  sulfur,  BTU/LB,  etc.   For  example: 

SELECT  mapname  ATTRIBUTE 

KEY?     ASH 

2  =  BOUNDED  SEARCH 

LOWER  BOUND  =  .01 

UPPER  BOUND  =  10  0 
GRID  the  attribute  ash 
PLC"  the  grid 

Repeat  -:hiH  procedure  for  each  attribute  recmired.   These  grids  car  be 
CONTOL"ed  and  SMOOTHed. 

SECTION  SIXTEEN:  PREPARING  THE  FINAL  MAPS 

For  screen  display: 

The  TEXT  command  creates  a  text  file  that  can  be  plotted  over  the 
maps  already  PLOTted  onto  the  screen.   The  AUTOLABEL  command 
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FIGURE   11.      Mine-striDoina   ratio  mao   aenerated  bv   the  QUINTIC  SPLINE 


inserts  labels  (values)  in  the  contour  lines  of  the  contour  map 
(line  map,  not  the  grid) . 

For  the  CALCOMP  or  similar  plotter: 

The  PENPLOT  command  produces  a  plot  file  of  a  map  at  a  user- 
specified  map  scale.   This  file  can  then  be  used  to  produce 
a  paper  or  mylar  map  on  a  plotter. 

SECTION  SEVENTEEN;  EXAMPLES  OF  ANALYSES 

FIGURES  12  through  14  show  various  ways  the  generated  maps  can  be 
used  in  MOSS  and  MAPS  to  analyze  data: 
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FIGURE    13 . 


Outcrop  map  generated  using   the   DEM    (Digital   Elevation 
Model)    for  the   surface  grid   and  substracting   the   structure 
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